Dynamical model for the quantum-to-classical crossover of shot noise 
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We use the open kicked rotator to model the chaotic scattering in a ballistic quantum dot cou- 
pled by two point contacts to electron reservoirs. By calculating the system-size-over-wave-length 
dependence of the shot noise power we study the crossover from wave to particle dynamics. Both 
a fully quantum mechanical and a semiclassical calculation are presented. We find numerically in 
both approaches that the noise power is reduced exponentially with the ratio of Ehrenfest time and 
dwell time, in agreement with analytical predictions. 
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I. INTRODUCTION 

Noise plays a uniquely informative role in connection 
with the particle-wave dualityi This has been appreci- 
ated for light since Einstein's theory of photon noise. 
Recent theoretical 2 * 3 ^*^ and experimental'' work has 
used electronic shot noise in quantum dots to explore 
the crossover from particle to wave dynamics. Particle 
dynamics is deterministic and noiseless, while wave dy- 
namics is stochastic and noisy£ 

The crossover is governed by the ratio of two time 
scales, one classical and one quantum. The classical time 
is the mean dwell time td of the electron in the quantum 
dot. The quantum time is the Ehrenfest time te, which is 
the time it takes a wave packet of minimal size to spread 
over the entire system. While td is independent of h, the 
time te increases oc ln(l/7i) for chaotic dynamics. An ex- 
ponential suppression oc exp(— te/td) of the shot noise 
power in the classical limit Ti — > (or equivalently, in the 
limit system-size-over-wave-length to infinity) was pre- 
dicted by Agam, Aleiner, and Larking A recent exper- 
iment by Oberholzer, Sukhorukov, and Schonenberger 7 
fits this exponential function. However, the accuracy 
and range of the experimental data is not sufficient to 
distinguish this prediction from competing theories (no- 
tably the rational function predicted by Sukhorukov 9 for 
short-range impurity scattering). 

Computer simulations would be an obvious way to test 
the theory in a controlled model (where one can be cer- 
tain that there is no weak impurity scattering to compli- 
cate the interpretation). However, the exceedingly slow 
(logarithmic) growth of te with the ratio of system size 
over wave length has so far prevented a numerical test. 
Motivated by a recent successful computer simulation of 
the Ehrenfest-time dependent excitation gap in the su- 
perconducting proximity effect^ we use the same model 
of the open kicked rotator to search for the Ehrenfest- 
time dependence of the shot noise. 

The reasoning behind this model is as follows. The 
physical system we seek to describe is a ballistic (clean) 
quantum dot in a two-dimensional electron gas, con- 



nected by two ballistic leads to electron reservoirs. While 
the phase space of this system is four-dimensional, it can 
be reduced to two dimensions on a Poincare surface of 
sectioniiiii^ The open kicked rotator^Sii 3 ^^ is a stro- 
boscope model with a two-dimensional phase space that 
is computationally more tractable, yet has the same phe- 
nomenology as open ballistic quantum dots. 

We study the model in two complementary ways. First 
we present a fully numerical, quantum mechanical solu- 
tion. Then we compare with a partially analytical, semi- 
classical solution, which is an implementation for this 
particular model of a general scheme presented recently 
by Silvestrov, Goorden, and one of the authors^ 



II. DESCRIPTION OF THE MODEL 

We give a description of the open kicked rotator, both 
in quantum mechanical and in classical terms. 



A. Closed system 

We begin with the closed system (without the leads). 
In this subsection we follow Refs. ITalT^ . The quantum 
kicked rotator has Hamiltonian 



H = 



KL 



h 2 d 2 



2I d0 2 T 



-cos 9 S s {t-kT ). (2.1) 



k= — oc 



The variable 9 € (0, 27r) is the angular coordinate of 
a particle moving along a circle (with moment of iner- 
tia Iq), kicked periodically at time intervals r (with a 
strength cx K cos 9) . To avoid a spurious breaking of 
time-reversal symmetry later on, when we open up the 
system, we represent the kicking by a symmetrized delta 
function: S s (t) = \8{t — e) + ^S(t + e), with infinitesi- 
mal e. The ratio %tq/2 , kIq = h e s represents the effective 
Planck constant, which governs the quantum-to-classical 
crossover. The stroboscopic time To is set to unity in 
most of the equations. 
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The stroboscopic time evolution of a wave function is 
given by the Floquet operator T = Tcxp(— i dt H/h), 
where T indicates time ordering of the exponential. For 
l//i e ff = M an even integer, T can be represented by an 
M x M unitary symmetric matrix. The angular coordi- 
nate and momentum eigenvalues are 8 m — 2irm/M and 
J m = n£, with m,i = 1,2,.. . M. We will use rescaled 
variables x = 0/2n and p = J/hM in the range (0, 1). 

The eigenvalues exp(— ie m ) of T define the quasi- 
energies e m £ (0,27r). The mean spacing 2-k/M of the 
quasi-energies plays the role of the mean level spacing 
5 in the quantum dot. In coordinate representation the 
matrix elements of T are given by 



mm/ • 

jj 1/2 2irimm r /M 

mm' — J'J ' 



(2.2a) 

e" '-, (2.2b) 

-i(MK/4ir) cos(27rm/M) in n \ 
mm' ^ , ^Z.ZCJ 

mm ,e-™ m2 ' M . (2.2d) 



The matrix product U^HU can be evaluated in closed 
form, resulting in the manifestly symmetric expression 



(U^UU) r 



M -i/2 e -i«/4 ex p[ i (7r/Af)(m' - mf 



(2.3) 

Classically, the stroboscopic time evolution of the 
kicked rotator is described by a map on the torus 
{x,p | modulo 1}. The map relates Xk+i,Pk+i at time 
k + 1 to Xk,Pk at time k: 



Xk+l 



Pk+i 



K 

Xk + Pk + -r- sm27ra;fe, 
4ir 



Pk 



— sin znxk 

47T 



sin 2ttx 



k+l 



(2.4a) 
(2.4b) 



The classical mechanics becomes fully chaotic for K > 7, 
with Lyapunov exponent A ~ hi(K/2). For smaller K the 
phase space is mixed, containing both regions of chaotic 
and of regular motion. We will restrict ourselves to the 
fully chaotic regime in this paper. 

For later use we give the monodromy matrix 
M(xk,Pk), which describes the stretching by the map 
of an infinitesimal displacement Sxk , Spk ■ 



Sxk+i 
Spk+i 



M(x k ,Pk)( S 6 X p k k 



From Eq. I|2.4|l one finds 

A-/Y \-( A ( Xk "> 1 

M{x k ,Pk)- ^ A (x fe )A(ar fe+1 ) - 1 A(x fe+1 ) 



K 

A(x) = 1 + — cos 2nx. 



(2.5) 

,(2.6a) 
(2.6b) 



B. Open system 

We now turn to a descriptio n of the open kicked ro- 
tator, following Refs. lldll5ll8L To model a pair of N- 
mode ballistic leads, we impose open boundary condi- 
tions in a subspace of Hilbert space represented by the 



indices in coordinate representation. The subscript 
n = 1,2, ...N labels the modes and the superscript 
a = 1,2 labels the leads. A 2N x M projection ma- 
trix P describes the coupling to the ballistic leads. Its 
elements are 



p — 

1 n.m. — 



1 if m = n £ {m 
otherwise. 



(«.) 



(2.7) 



The matrices P and T together determine the quasi- 
energy dependent scattering matrix 



S(e) = P[e- l£ 
Using PP T = 1, Eq. 

PAP T - 1 



5' 



P T P)Y 1 TP T . (2.8) 
can be cast in the form 
1 + e ie T 



PAP T + 1 ' 



A 



1 - e ie T 



-A\ 



(2.9) 



which is manifestly unitary. The symmetry of J- ensures 
that S is also symmetric, as it should be in the presence 
of time-reversal symmetry. 

By grouping together the N indices belonging to the 
same lead, the 2N x 27V matrix S can be decomposed 
into 4 sub-blocks containing the N x N transmission and 
reflection matrices, 



S 



r t 
t' r' 



The Fano factor F follows fromi^ 

Tr«t(i _ tt t) 



F = 



Trttt 



(2.10) 



(2.11) 



This concludes the description of the stroboscopic 
model studied in this paper. For completeness, we briefly 
mention how to extend the model to include a tunnel bar- 
rier in the leads. 

To this end we replace Eq. 1)2. 8fl by 

S(e) = -(1-KK T ) 1 / 2 
+ K[e- l£ - T(l - K^Kfl^TK 1 \ (2.12) 



The N x M coupling matrix K has elements 



Kn 



yTV^ if m = n £ {m' r . 
otherwise, 



(2.13) 



with r„ £ (0, 1) the tunnel probability in mode n. Ballis- 
tic leads correspond to T n = 1 for all n. The scattering 
matrix l|2.12|l can equivalently be written in the form 
used conventionally in quantum chaotic scattering: 20 ' 21 



S(e) = -1 + 2W(A~ 1 + W rT W)' 1 W T , 



(2.14) 



with W = K(l + Vl - K^K)- 1 and A defined in Eq. 
(EU). 
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FIG. 1: Dependence of the Fano factor F on the dimensionality of the Hilbert space M — l/h e s, at fixed dwell time td = M/2N 
and kicking strength K. The data points follow from the quantum mechanical simulation in the open kicked rotator. The solid 
line at F = \ is the M-independent result of random-matrix theory. The dashed lines are the semiclassical calculation using 
the theory of Ref. |f| There are no fit parameters in the comparison between theory and simulation. 



III. QUANTUM MECHANICAL CALCULATION 

To calculate the transmission matrix from Eq. I|2.8[l we 
need to determine an NxN submatrix of the inverse of an 
M x M matrix. The ratio M/2N — rrj is the mean dwell 
time in the system in units of the kicking time tq. This 
should be a large number, to avoid spurious effects from 
the stroboscopic description. For large M/N we have 
found it efficient to do the partial inversion by iteration. 
Each step of the iteration requires a multiplication by 
which can be done efficiently with the help of the fast- 
Fourier-transform algorithmi22i22i We made sure that the 
iteration was fully converged (error estimate 0.1%). In 
comparison with a direct matrix inversion, the iterative 
calculation is much quicker: the time required scales oc 
M 2 In M rather than oc M 3 . 

To study the quantum-to-classical crossover we re- 
duce the quantum parameter h e s — 1/M by two or- 
ders of magnitude at fixed classical parameters td = 
M/2N = 5, 10,30 and K = 7, 14,21. (These three val- 
ues of K correspond, respectively, to Lyapunov expo- 
nents A = 1.3,1.9,2.4.) The left edge of the leads is at 
m/M = 0.1 and m/M = 0.8. Ensemble averages are 
taken by sampling 10 random values of the quasi-energy 
e £ (0, 2tt). We are interested in the semiclassical, large- 
TV regime (typically N > 10). The average transmission 
N^ 1 (Trtt' i ) w 1/2 is then insensitive to the value of h e g, 
since quantum corrections are of order 1/N and there- 
fore relatively smalliSi The Fano factor l|2.11|l . however, 
is seen to depend strongly on h e s, as shown in Fig. ^ 
The line through the data points follows from the semi- 
classical theory of Ref.|| as explained in the next section. 

In Fig. we have plotted the numerical data on a 
double-logarithmic scale, to demonstrate that the sup- 
pression of shot noise observed in the simulation is in- 
deed governed by the Ehrcnfcst time te ■ The functional 
dependence predicted for N > \[M is£ 
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FIG. 2: Demonstration of the logarithmic scaling of the Fano 
factor F with the parameter N 2 /M = M/(2t d ) 2 . The data 
points follow from the quantum mechanical simulation and 
the lines are the analytical prediction 13. II . with c a fit pa- 
rameter. The slope A -1 = l/ln(if/2) of each line is not a fit 
parameter. 



F 



-te/td 



X- 1 \n(N 2 /M) + c, (3.1) 



with c a if-dependent coefficient of order unity. As 
shown in Fig. |2J the data follows quite nicely the log- 
arithmic scaling with N 2 /M = M/(2te>) 2 predicted by 
Eq. H3.1fl . This corresponds to a scaling with w 2 / LAp in 
a two-dimensional quantum dot (with A p the Fermi wave 
length and w and L the width of the point contacts and 
of the dot, respectively.) We note that the same para- 
metric scaling governs the quantum-to-classical crossover 
in the superconducting proximity effect ^2*24 
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FIG. 3: (Color online.) Upper panels: phase space portrait of lead 1, for td = 10 and different values of K. Each point 
represents an initial condition for the classical map 12.41 . that is either transmitted through lead 2 (black/red) or reflected 
back through lead 1 (gray/green). Only initial conditions with dwell times < 3 are shown for clarity. Lower panels: histogram 
of the area distribution of the transmission and reflection bands, calculated from the corresponding phase space portrait in the 
upper panel. Areas greater than h e g = 1/M correspond to noiseless scattering channels. 



IV. SEMICLASSICAL CALCULATION 



To describe the data from our quantum mechanical 
simulation we use the semiclassical approach of Ref. |3 
To that end we first identify which points in the x,p 
phase space of lead 1 are transmitted to lead 2 and which 
are reflected back to lead 1. By iteration of the classical 
map H2.4f) we arrive at phase space portraits as shown 
in Fig. |31 (top panels). Points of different color (or gray 
scale) identify the initial conditions that are transmitted 
or reflected. 

The transmitted and reflected points group together 
in nearly parallel, narrow bands. Each transmission or 
reflection band (labeled by an index j) supports noiseless 
scattering channels provided its area Aj in phase space 
is greater than h B s = 1/M. The total number iVo of 



noiseless scattering channels is estimated by 

3 

with 6[x) = if x < and 6(x) = 1 if x > 0. In the 
classical limit M — ► oo one has Nq = N, so all channels 
are noiseless and the Fano factor vanishes £ 

As argued in Ref.0, the contribution to the Fano factor 
from the N — No noisy channels can be estimated as 
1/4N per channel. In the quantum limit No = one 
then has the result F = 1/4 of random-matrix theory^ 
The prediction for the quantum-to-classical crossover of 
the Fano factor is 

3 

M f 1/M 

= AN Jo A P^ dA > ( 42 ) 






0.1 



1 
0.15 0.1 





0.15 0.1 



0.15 



FIG. 4: (Color online.) Contour plots of the Husimi function ffHty in lead 1 for M = 2400, r D = 10, and K = 7, 14, 21. The 
outer contour is at the value 0.15, inner contours increase with increments of 0.1. Yellow regions are the classical transmission 
bands with area > 1/M, extracted from Fig. |3] 



with band density p(A) = Y], 5 (A — Aj). The quantum 

limit F = 1/4 follows from the total area f Q Ap(A) dA = 
N/M. 

We have approximated the areas of the bands from the 
monodromy matrix (|2.6[) . as detailed in the Appendix. 
The lower panels of Fig. [3] show the band density in the 
form of a histogram. The solid curves in Fig. ^ give the 
resulting Fano factor, according to Eq. 1)4.2(1 . 

V. SCATTERING STATES IN THE LEAD 

To investigate further the correspondence between the 
quantum mechanical and semiclassical descriptions we 
compare the quantum mechanical eigenstates \Ui) of t'H' 
with the classical transmission bands. 

Phase space portraits of eigenstates are given by the 
Husimi function 

Hi(m x ,m p ) = \(Ui\m x ,m p )\ 2 . (5.1) 

The state \m x ,rn p ) is a Gaussian wave packet centered 
at x = m x /M, p — rrip/M. In position representation it 
reads 

oc 

{m\m x ,m p ) oc y, & e p ' . 

k— — co 

The summation over k ensures periodicity in m. 

The transmission bands typically support several 
modes, thus the eigenvalues Ti are nearly degenerate at 
unity. We choose the group of eigenstates with Ti > 
0.9995 and plot the Husimi function for the projection 
onto the subspace spanned by these eigenstates: 

li(m x ,m p ) = 2_j T^i{ m x,m p ). (5.3) 

T;> 0.9995 

As shown in Fig. 0] this quantum mechanical function 
indeed corresponds to a phase-space portrait of the clas- 
sical transmission bands with area > 1/M. 



VI. CONCLUSION 

We have presented compelling numerical evidence for 
the validity of the theory of the Ehrenfest-time dependent 
suppression of shot noise in a ballistic chaotic system m& 
The key prediction^ of an exponential suppression of the 
noise power with the ratio te /td of Ehrenfest time and 
dwell time is observed over two orders of magnitude in the 
simulation. We have also tested the semiclassical theory 
proposed recently^ and find that it describes the fully 
quantum mechanical data quite well. It would be of in- 
terest to extend the simulations to mixed chaotic/regular 
dynamics and to systems which exhibit localization. 
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APPENDIX: CALCULATION OF THE BAND 
AREA DISTRIBUTION 

We approximate the bands in Fig. 01 by straight and 
narrow strips in the shape of a parallelogram, disregard- 
ing any curvature. This is a good approximation in par- 
ticular for the narrowest bands, which are the ones that 
determine the shot noise. Each band is characterized by 
a mean dwell time n (in units of To). We disregard any 
variations in the dwell time within a given band, assum- 
ing that the entire band exits through one of the two leads 
after n iterations. (We have found numerically that this 
is true with rare exceptions.) 
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FIG. 5: Phase space of a lead (width w) showing two areas 
(in the shape of a parallelogram) that are mapped onto each 
other after n iterations. They have the same base B, so the 
same area, but their tilt angle a is different. 

The case of a reflection band is shown in Fig. The 
initial and final parallelograms have the same height, set 
by the width w = N/M of the lead. Since the map is area 
preserving, the base B of the two parallelograms should 



be the same as well. To calculate the band area A = Bw 
we assume that the monodromy matrix M(xk,Pk) does 
not vary appreciably within the band at each iteration 
k = 1,2, .. .n. An initial vector e*j within the parallelo- 
gram is then mapped after n iterations onto a final vector 
ef given by 

e) = Me h M = M{x n ,p n ) ■ ■ ■ M(x 2 ,P2)M(xi,px), 

with xi,pi inside the initial parallelogram. 

We apply Eq. I|A.1|) to the vectors that form the sides 
of the initial and final parallelograms. The base vector 
e - t = Bp is mapped onto the vector ej = ±w(i+ptana), 
with a the tilt angle of the final parallelogram. It follows 
that B\M xp \ — w, hence 

A = w 2 /\M xp \. (A.2) 

We obtain the Fano factor F by a Monte Carlo pro- 
cedure. An initial point x\,p\ is chosen randomly in 
lead 1 and iterated until it exits through one of the two 
leads. The product M. of monodromy matrices starting 
from that point gives the area A of the band to which it 
belongs, according to Eq. (|A.2|) . The fraction of points 

with A < 1/M then equals w' 1 J Q 1/M Ap(A) dA = 4F, 
according to Eq. (|4.2[) . 

To assess the accuracy of this procedure, we repeat the 
calculation of the Fano factor with initial points chosen 
randomly in lead 2 (instead of lead 1). The difference is 
about 5%. The dashed lines in Fig. Q are the average of 
these two results. 
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